Relative production rates of °He, ^Be, C in astrophysical environments 
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We assume an environment of neutrons and a-particles of given density and temperature where 
nuclear syntheses into ^He, ®Be and ^^C are possible. We investigate the resulting relative abundance 
as a function of density and temperature. When the relative abundance of a-particles Ya is between 
0.2 and 0.9, or larger than 0.9, the largest production is ^Be or ^'^C, respectively. When Ya < 0.2 
^He is mostly frequently produced for temperatures above about 2 GK whereas the ®Be production 
dominates at smaller temperatures. 

PACS numbers: 25.10.+S, 25.40.Lw, 26.30.Hj 



Introduction. Formation of heavy elements must 
overcome the problem that all nuclear isotopes with mass 
numbers 5 and 8 are unstable [l[ ■ As a consequence, once 
the hydrogen fuel in a star is exhausted, the production 
of energy by formation of ^He stops and the temperature 
drops. The subsequent gravitational collapse increases 
the temperature, and the red giant phase, where helium 
is now the source of energy, begins. Due to the lack of 
neutrons in the core of the helium burning red giants, the 
A=5 and A=8 instability gaps have to be bridged by the 
triple-alpha reaction a + a + a —^ ^^C -I- 7, which is the 
most relevant one in stars at the helium burning stage. 

Nevertheless, at conditions of high alpha and neutron 



densities, the reactions a 



^He -|- 7 and a 



a + n ^ ■^Be + 7 are also possible, and they will play a 
role in bridging the A = 5,8 gaps. As described in [ll-Q, 
such a scenario can appear at the early Big Bang stages 
or in the nucleosynthesis related to the type II supernova 
shock front. In both cases, temperatures are estimated 
to be of about 7 to 10 GK (1 GK=10S K). 

The dominating process is in any case the electromag- 
netic radiative recombination of the three particles from 
continuum to bound state, except at very high densities 
where four-body recombination can compete favorably 
[J|. However, the three-body processes producing these 
nuclei differ from each other. In particular, the produc- 
tion of ^He and ^Be is dominated by dipole transitions, 
while the production of ^^C is of quadrupole character. 

The three-body processes described above are usually 
treated as two-step processes where the unstable isotopes 
^He and *Be first are created and second, before decay- 
ing, react with another neutron or a-particle [5|-l7|. This 
approximation of two independent sequential processes 
does not provide an accurate description since the life- 
times of the intermediate configurations are comparable 
to, or shorter than, the reaction time of the last step of 
the process Q. This implies that the processes proceed 
through genuine three-body reactions. 

The approximations employed to overcome this diffi- 
culty are, almost always, extensions of the bare two-step 
model where the finite lifetime, or width, of the reso- 



nances is incorporated into the description [9|. The fact 
that this can be tremendously wrong was discussed re- 
cently in a full three-body calculation of the triple a- 
process at very low energies. In [lOj an increase in the 
triple alpha reaction rate by about 20 orders of magni- 
tude around 10^ K compared with the rate in [ll| was 
found. In [12] it is shown that such an increase is incom- 
patible with observations of extended red giant branches 
and He burning stars in old stellar systems. The three- 
body character is in general crucial for resonance struc- 
ture calculations 13, 14] and the related dynamic evolu- 
tion describing the decay mechanism [14l - ll6| . 

The correct three-body calculations are technically dif- 
ficult because they involve the continuum dynamics of 
three particles interacting via a mixture of short and 
long-range forces. Only recently it has become feasible to 
perform general three-body computations with the cor- 
rect boundary conditions at both small and large dis- 
tances, and where no assumption is made about the cap- 
ture mechanism (sequential or direct). 

The purpose of the present letter is to investigate 
the only possible three-particle reactions at conditions 
of high alpha and neutron density. The novelty of this 
work is that we are doing it using a three-body method 
that is not making any assumption about the reaction 
mechanism. In this way, we shall investigate the rela- 
tive abundances of the three nuclei, ^He, ^Be and ^^C, 
created from neutrons and a-particles at given densities 
and temperatures. The conditions under which the pro- 
duction of each nucleus is the dominant process will be 
established. 

Formulation. We consider radiative capture of three 
particles, (abc), into a bound nucleus A of binding energy 
Ba, i.e. a + b-\-c^A + j. The particles are a-particles 
and neutrons in various combinations. The correspond- 
ing reaction rate Rabc{E) can be related to the inverse 
process of photodissociation. Eq.(20) in jlTf gives the ra- 
tio between the "energy averaged" reaction rates for the 
processes a + b + c^^A + d. Making use of (32) in [171 
and specializing to the case in which d is a photon, one 



can identiiy the relation: 
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where E — E^ + Ba is the initial three-body kinetic en- 
ergy, E~f is the photon energy, a^{E^) is the photo disso- 
ciation cross section of the A nucleus, c is the velocity of 
light, gi is the degeneracy of particle i = a,b, c, A, and fj,x 
and Hy are the reduced masses of the systems related to 
the Jacobi coordinates, (a;, y), for the three-body system 

The production rate P for the capture reaction is 
obtained after averaging Rabc{E) using the ]VIaxwell- 
Boltzmann distribution as weighting function, and mul- 
tiplying by the density rii of particles o, b, and c. This 
density is usually written as n, = pNAXi/Ai, where 
p is the density of the environment, Na is the Avo- 
gadros number, Ai is the mass number of particle z, and 
Xi = NiMi/{NaMa + NbMb + NcMc) is the mass abun- 
dance of nucleus i expressed by the number of particles 
Ni and their masses Mi (see Eqs.(l) and (3) in ([i4|)- It 
is also possible to use the relative abundance defined by 
Y, = Nj{Na + Nb + N,). 

In this way the final expression for the production rate 
P depends then on both, temperature (T) and mass den- 
sity (p) of the environment, which can vary substantially 
in different scenarios, and it takes the form fl7| : 



Pabc{p,T) = Uamn 
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The photodissociation cross section for the inverse pro- 
cess A + ^ ^ a-\-h + c can be expanded into electric and 
magnetic multipoles. In particular, the electric multipole 
contribution of order A has the form |19i: 
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where the strength function B is 
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where Jq, J and Mq, M are the total angular momenta 
and their projections of the initial and final states, and 
all the other needed quantum numbers are collected into 
no and n. The electric multipole operator is given by: 



o^^ = J2z^\n-R\^Yx^,iny 



(5) 



where i runs over the three particles, and where we ne- 
glect contributions from intrinsic transitions within each 
of the three constituents ,20,1 . 



Test of the method. The crucial strength function, B, 
is computed by genuine three-body calculations of both 
the bound final state and the continuum initial states. 
We have used the hyperspherical adiabatic expansion 
method described in |18| . The n-n, a-n, and a-a in- 
teractions are the ones given in [ij, [21| . The basic pro- 
cedure is computation of three-body states of given an- 
gular momentum and parity confined by box boundary 
conditions [23|- In this way the continuum spectrum is 
discretized. The strength functions are then obtained for 
each discrete continuum state according to Eq.®. The 
distribution dB/dE is built by use of the finite energy in- 
terval approximation, where the energy range is divided 
into bins, and all the discrete values of B falling into a 
given bin are added. Afterwards the points are connected 
by spline operations and the expressions ([3]) and ^ are 
computed. 
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FIG. 1: (a) Photodissociation cross section for 'Be. The total 
computed cross section is given by the solid line. The contri- 
butions from transitions to the 1/2"*" , 3/2^ , and 5/2^ contin- 
uum three-body states are given by the dashed, dotted, and 
dot-dashed curves, respectively. The experimental data are 
from [g. (b) Comparison between different computed dipole 
strength functions for ^He and the experimental data (shaded 
area) from 23]. 

In the upper part of Fig [1] we show the computed pho- 
todissociation cross section of ^Be. The experimental 
data in the figure are from [y]. The dashed, dotted, and 
dot-dashed curves give the contribution to the transition 



strength in Eq.Q, and therefore to the total cross sec- 
tion (sohd hne), from the a + a + n continuum states with 
total spin and parity 1/2+, 3/2+, and 5/2+, respectively 
(the ground state in ^Be has spin and parity 3/2~). As 
seen in the figure, the agreement with the experiment is 
reasonably good, except for energies above 2 MeV, where 
the computed cross section clearly overestimates the ex- 
perimental data. This is because the interaction used 
in the three-body calculation provides a 3/2+ resonance 
twice wider than the experimental one. The narrow peak 
at 0.9 MeV is a known 5/2~ resonance which can only be 
populated through the Ml transition which is neglected 
in this calculation. 

For ^He, our computed dB/dE (thick solid line in the 
figure) agrees reasonably well with the previous calcu- 
lations in [2J, [2^. The one in [26| gives clearly larger 
values, mainly for energies beyond 1 MeV. In any case, 
none of the calculations fits well the experimental data 
(thin solid line) for all the energies despite the large ex- 
perimental error bars. All the calculations show a sharp 
enhancement in the low energy region which is not seen 
in the experiment. Further experimental data would help 
to clarify this issue. In similar measurements with ^^Li, 
the peak in the ex per imental strength function moved to 
lower energies_ 23, [2g| , and thus closer to our theoretical 
predictions [29|, after improvements in the detection of 
the low energy neutrons. 

Density dependence. The density dependence of the 
production rates in Eq. ^ is very simple for a given tem- 
perature. The basic reaction rate for only three particles 
has to be multiplied by the number of each species of 
particles. For a given total density p found by adding 
neutron and a-particle densities we can express the den- 
sity dependence as p'^X^AT^"", where Xn = 1 — Xa and 
Xa are the mass fraction of neutrons and a-particles. 
Then n = 1,2,3 correspond to production of ^He, ^Be 
and ^^C, respectively. 

When no a-particles are present, Xa = Yq = 0, the 
production rates are all zero. When only a-particles are 
present, X^ — Ya — 1, only ^^C can be produced. The 
density dependence for production of ^He and ^Be are 
each others reflection in Xa = 1/2 but pushed towards 
smaller values as function of Yq. The production of ^^C 
increases monotonically as function of Xa and Ya ■ 

Temperature dependence. The temperature depen- 
dence is obtained by folding the calculated energy de- 
pendence with the Boltzmann distribution as seen for- 
mally in Eq.(I2]) and shown numerically in Figj^l We 
show results from about 0.1 GK to 10 GK. The lower 
limit is where we run out of accuracy due to the present 
discretization method, and the upper limit allows inclu- 
sion of the relevant processes. For completeness we dis- 
cuss the behavior in the full interval. The three-body 
processes producing the three nuclei, ^He, ^Be, ^^C, dif- 
fer from each other. Two neutrons and one a-particle 
produce the 0+ ground state of ^He by electromagnetic 
recombination where one photon emerges with total an- 
gular momentum and parity J'^ = 1~,2+. The dipole 
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FIG. 2: The rates as functions of temperature for producing 
He, Be and C in their respective ground states from var- 
ious continuum states. From 1~ and 2+ to 0+ for ^He (solid 
curves), from 1/2+ , 3/2+, and 5/2+ to 3/2" for '^Be (dashed 
curves), from 2+ to the 0+ ground state and from both 0+ 
and 2+ to 2+ excited bound state for ^^C (dotted curves). 



transition is strongly preferred in spite of the low-lying 
2+ three-body resonance in ^He [4|.|23|. 

One neutron and two a-particles in continuum states of 
J^ = 1/2+, 3/2+, 5/2+ combine directly by dipole tran- 
sitions into the 3/2~ ground state of ^Be. Resonances 
enhance the transitions in various energy ranges. At low 
temperature transitions from the 1/2+ state dominate by 
several orders of magnitude but at T sa 5 GK the 5/2+ 
transition has almost reached the same value, see Figl2] 
This behavior is related to the corresponding resonances 
at excitation energies of 1.68 MeV and 3.05 MeV, respec- 
tively. The 3/2+ resonance is located at a much higher 
energy and the corresponding transition always remains 
smaller. Furthermore, higher multipolcs give insignifi- 
cant contributions at these temperatures |^. 

Three a-particles can combine into the ground state 
of ^^C through several paths. At low temperature the 
quadrupole transitions from 0+ dominates due to the 
Hoyle resonance, but the quadrupole transition from the 
2+ continuum to the ground state increases and begin 
to dominate at T « 2.5 GK, see Figl^l The precise 
value depends on the energy of the lowest 2+ resonance 
which here was assumed to be 1.38 MeV above the three- 
body threshold |14| . The remaining quadrupole transi- 
tion from 2+ to 2+ is smaller for all temperatures. 

All these transitions are compared in Figl^l The pro- 
duction rate for ^Be is by far the largest at low tem- 
peratures but matched for ^He above temperatures of 
T « 8 GK. The rate for production ^^C is much smaller 
in this temperature range. 

In Fig[3]the total reaction rates in Figl5]are compared 
with previous calculations. The computed rates for ^He, 
^Be, and ^^C are given by the dashed, sohd, and dotted 
curves, respectively. In [3| the reaction rate for produc- 
tion of ^He is given by assuming a sequential capture 
process through the 3/2~ resonance in ^He. This reac- 
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FIG. 3: Comparison between the total computed rates for 
production of He (dashed line), Be, (solid line), and ^ C 
(dotted line) and the ones calculated in [7| for ®He, in [g, |3Q] 
for ^Be, and in 0,113 for ^^c. 



tion rate is about three orders of magnitude smaller than 
our calculation, where the sequential capture assumption 
has not been made. In [3] they give the same reaction 
rate but adding the contribution of a dineutron capture. 
This is shown by the closed triangles in the figure, which 
now an order of magnitude above our calculation. 

For production of ^Be, our reaction rate agrees well 
with the ones in Refs.|3(]{ (open squares) and Q (closed 
squares) for temperatures above 1 GK. Below this tem- 
perature our reaction rate is clearly bigger. For ^^C our 
computed rate agrees very well with the results given in 
[111 |30| , where a sequential capture through the narrow 
0+ resonance in ^Be is assumed. The agreement with [11] 
(closed circles) is particularly good. The deviations from 
[30] (open circles) arise because we include the contribu- 
tion from the 2"*" — >■ 0^ transitions which is responsible 
for the enhancement of the tail for temperatures above 3 
GK (see the dotted curves in Figl2]). The good agreement 
with 11] does not imp ly that we disagree with the spec- 
tacular results in [lO| where huge enhancements of the 
reaction rates are reported for energies around 0.01 GK. 
We only discuss temperatures down to 0.1 GK where the 
present method is applicable. However, we see no indi- 
cation of such rate increase at very low temperature. 

Density-temperature dependence. The complicated 
temperature dependence in Figl2] can now be combined 
with the simple density dependence, nanbUc, in Eq.(I2]). 
After adding up the different contributions for each of the 
three nuclei we find the creation probabilities as functions 
of temperature for three different relative fractions of the 
three nuclei. Creation of ^^C only becomes competitive 
for any temperature when Yq is close to unity, i.e. when 
essentially only a-particles are present. 

The slightly increasing ^He curves cross the slightly de- 
creasing ^Be curves while the ^^G curves stay below for 
all temperatures below 10 GK. Thus unless the relative 
a-neutron abundance is extreme we get dominance of ^Be 
at small temperature and dominance of ^He at large tem- 



perature. However, now it is important precisely which 
number of neutrons and a-particles are available for the 
recombination process. When Yq — 0.1 the dominance 
changes from ^Be to ^He at T « 2.5 GK. For Ya values 
larger than 0.3 ^Bc dominates for all temperatures. 
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FIG. 4: The phase diagram for producing ^He, ^Be and ^^C in 
the ya-temperature parameter space. The curves correspond 
to a constant ratio of production rates of two nuclei. 



Relative production. The complete picture of the 
density-temperature dependence of the three-body reac- 
tion rates for creating the three nuclei is shown in FigU) 
All these nuclei can be destroyed in the same environment 
but we consider here only the individual three-body re- 
actions creating them. Transfer reactions from ^He and 
^Be using the same basic ingredients of neutrons and al- 
phas can also produce ^Be and ^^G. However, we leave 
these four-body processes for future work, they are be- 
yond the scope of the present investigation. For Ya less 
than about 0.1 and temperatures above T « 1 — 4 GK the 
nucleus ^He is produced more than twice as often as ^Be. 
As Ya increases the relative ^Be production increases and 
becomes dominant for all temperatures when Ya exceeds 
0.2. 

By further increase of Ya the relative creation rate of 
^^C increases. At Ya « 0.82 the ^Be production is still 
dominating when -'^^C is created with half the rate of ®He. 
At Ya « 0.86 the production rates for ^^C and ^He are 
equal, but the production of ^Be still dominates. Only 
when Ya is larger than about 0.99, where very few neu- 
trons are present, the production rate of "'^^C exceeds the 
other rates. These relative rates are very crudely inde- 
pendent of temperature except for very low Ya values. 
Except for the small Ya results, similar overall conclu- 
sions were obtained in previous investigations (j, [a, [3l[ . 

Summary and conclusions. We have compared calcu- 
lated results of production rates of three-cluster nuclei, 
^He, ^Be, ^^C, consisting of neutrons and a-particles. 
The density and temperature dependence are investi- 
gated in a three-body model where the necessary con- 
tinuum structures are treated by discretization in a large 
box. The two-body interactions are realistic and repro- 
ducing low-energy scattering properties. The position of 



the lowest 2+ resonance in ^^C is not precisely known but 
the results are rather insensitive to this energy. In gen- 
eral the results are very robust against variations within 
realistic limits of the interactions. 

The ^^C production only dominates when the num- 
ber of a-particles is 100 times larger than the number 
of neutrons. At lower relative a-particle density the ^Be 
production is most frequent until ®He begins to domi- 
nate when the number of a-particles is about 1/4 of the 
number of neutrons. An exception is low a-particle den- 
sity and temperatures below about 2 GK where the ^Be 
rate is largest. Thus the actual temperature is important 
when only relatively few a-particles are present. 



The full implications of the present results require de- 
tailed numerical investigations of the sequences of pro- 
cesses leading to formation of heavier nuclei. The initial 
conditions in these related chains of reactions should be 
our relative production rates. Whether the outcome is 
consistent with observations of nuclear abundances re- 
mains to be seen. 
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